rm(list = ls(all = TRUE)) #clear workspace

library(foreign)
library(survey)
options(survey.lonely.psu = "adjust")

setwd("C:/Users/martin.romero/Dropbox/ENVIPE/2012/BD_ENVIPE12")
#setwd("~/Dropbox/ENVIPE/2011/Base_de_datos")

################################################
#Pregunta al entrevistado
################################################

data <- read.dbf("tper_vic.dbf")

#Ponderador y diseno
data$FAC_ELE <- as.numeric(as.character(data$FAC_ELE))
design.per <- svydesign(id=~UPM,strata=~EST, weights=~FAC_ELE, data=data)

#Lo secuestraron para exigir dinero o bienes:
# 1 Si
# 2 No
# 8 No responde
# 9 No sabe
table(data$AP7_3_8) #78 casos
foo <- svymean(~AP7_3_8, design=design.per, na.rm=TRUE) #0.12%
print(foo)
confint(foo)
foo <- svytotal(~AP7_3_8, design=design.per, na.rm=TRUE) #93,262
print(foo)
confint(foo)

#Me puede decir cuantas veces?
table(as.numeric(as.character(data$AP7_4_8)))
sum(as.numeric(as.character(data$AP7_4_8)), na.rm=TRUE) #81
foo <- svytotal(~as.numeric(as.character(AP7_4_8)), design=design.per, na.rm=TRUE) #98,826
print(foo)
confint(foo)

################################################
#Modulo de victimizacion
################################################

modulo <- read.dbf("tmod_vic.dbf")

#6.2 �Me podr�a decir si sus secuestradores.  					
#1  exigieron rescate a familiares por su liberaci�n?
#2	"lo obligaron a retirar dinero de un cajero, entregar joyas, celular u otras cosas?"
#9	No sabe / no responde
#b	blanco
table(modulo$BP6_2)
round(prop.table(table(modulo$BP6_2)),3)*100

#6.1 �Cu�nto tiempo lo tuvieron secuestrado?  					
#1  Menos de 24 horas
#2	De 1 a 3 d�as
#3	De 4 a 10 d�as
#4	De 11 a 29 d�as
#5	De 1 a 3 meses
#9	No sabe / no responde
#b	blanco
table(modulo$BP6_1)
round(prop.table(table(modulo$BP6_1)),3)*100



